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1. Introduction 


Among the vast scope of matrix ensembles a distinguished role together with the in¬ 
tegrals over a single matrix is played by two matrix models - the ensemble of two matrices, 
usually with the simplest possible interaction between them []lj. Being still simple inte¬ 
grate systems (like in the one matrix case their partition function is a tau-function of 
Toda lattice hierarchy [Q]), these models possess already a richer mathematical structure 
than one matrix models and thus give rise to more applications. The two matrix model was 
proposed and studied in [Q as an important solvable example of a new class of statistical- 
mechanical systems: spins with nearest neighbor interaction on planar graphs (Ising spins 
in this case). Its multi-critical generalization, in the spirit of the one matrix multi-critical 
points |Q], leads to the complete picture of (p, g)-critical points in two-dimensional gravity 
It appears also in the context of two-dimensional Laplacian growth i,l, 0 
demonstrating some hidden parallels between all these problems. 

Matrix models have been the first physical example when the partition functions 
were directly related to tau-functions of integrable systems [!2[|. The relation between the 
(quantum) partition functions and tau-functions of (classical) integrable systems is still 
rather intriguing, it appears to be much more universal than one could expect before. 
The same sort of integrable structure like in the matrix models and/or the topological 
string theories has been found jTlJ in the context of Seiberg-Witten theories | p~2 or J\f — 2 
supersymmetric gauge theories in four dimensions. 

The similarity between matrix models and supersymmetric gauge theories based on 
the similarity of their integrable structures was noticed long ago ]I3|. However, in multi¬ 
dimensional supersymmetric gauge theories, apart from a recent example [[Uj, one mostly 
observed the quasiclassical limit of integrable hierarchies (see ]I5| for details and refer¬ 
ences). It means that the prepotentials of gauge theories are described rather in terms of 
the quasiclassical tau-functions or tau-functions of universal Whitham hierarchy |I6[], than 
the tau-functions of dispersionfull hierarchies. 

The overwhelming part of the work on the matrix models in the planar (large N ) 
limit concerned the so called one cut case, when the eigenvalues form a single support dis¬ 
tribution, though a few interesting papers on the multi-support distributions were written 
in the past, especially jT7j] (see also |0|), where the relation with the hyperelliptic curves 
was revealed. All these papers were devoted to the one matrix model. 


Recently, a new interest to the multi-cut solutions was born, due to the papers ||T9 
where the effective superpotentials in the J\f — 1 supersymmetric gauge theories were 
related to the matrix models. The multi-cut solution corresponds there to the breaking 
of the gauge group into a few subgroups. It was also proposed in [|II| to “fill” by the 


1 A similar picture arises in the generalized one matrix model in external field [7], and it is not 
surprising that the multi-support solutions in these models are also related. 
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eigenvalues not only the minima of the matrix potential, but also the maxima, the situation 
missed in the matrix model literature, due to the obvious absence of stability of such 
configurations. However, from the mathematical point of view, especially from the point 
of view of the analytic curve description such filling appears to be admissible and even 
having some nice physical applications. These aspects of the multi-cut solution were further 
developed in recent papers |T| , |[ZT]| , |[Z2] . 

In this paper we are going to study the multi-support solutions in the two matrix 
model. As in the one matrix model case, these solutions can be formulated in terms of 
geometry of the underlying complex curve endowed with a generating differential and a 
quasiclassical tau-function, of the type proposed in |16| . We will write down explicitly the 
equation for the complex curve in the two matrix model and define the partition function 
in the planar limit as a quasiclassical tau-function. We will also study an important degen¬ 
erate case of real potentials and demonstrate the consistency of the tau-function approach 
with the planar graph expansion (in terms of a multi-phase Ising system on the graphs) 
for the multi-support solution in the two matrix model. Our method of construction of 
the algebraic curve of the two matrix model is powerful enough to be generalizable to the 
more general multi-matrix models with the so called tree-like interactions of matrices. 

The two matrix model with the simplest interaction between matrices is known in the 
literature in two, superficially different, forms. The partition function of the normal matrix 
model (NMM)i of two commuting complex conjugated N x N matrices $, <fd: [$, <fd] = 0 
is defined as follows: 


Z N [t,t] = J V $ V& e -Ti 'K$ t +2ReT 1 'lV($)_ 


( 1 . 1 ) 


where the harmonic part of the potential H($, <f>f) = —+ 2Re W(5>) is parameterized 
as W(<f>) = Yhk=i ■ Going to the eigenvalues $ = diag(zi,..., zn) we obtain!: 

r N 

Z i * MM) [t,i]= I JJ (d 2 3 z m e -*mZm+2ReW(z m )^ |A(z)| 2 , (1.2) 

J m= 1 V 

the normalization by the unitary group volume Vjjr jyj being hidden into the definition of 
integration measure in (0>. The last integral has a natural interpretation in terms of 
the partition function of Coulomb gas of particles with coordinates Zj, confined by the 
potential in the exponent. 

Almost the same eigenvalue integral can be presented as the partition function of the 
model of two hermitian non-commuting matrices X , Y (H2MM) 


4f 2MM) [t,(]= / VXVY e --nXY+T,W(X)+TtW(Y) 


(1.3) 


2 which is a particular case of the models defined in [23] 

3 after taking into account the Jacobian of the angular part of commuting matrices 
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where W (X) = J^k=i tkX k and W ( Y ) = Yhk=i ^kY k or, in terms of eigenvalues (using, as 
usual, the Harish-Chandra-Itzykson-Zuber (HCIZ) formula): 




(H2MM) 

N 


r N 

[t,t}= / n (dx m dy m e -^y m +w(x m )+w(y m )^ A(x)A( , y) (L4) 

^ m=1 


Now the Jacobian gives X(.x) J A(y) 2 in contrast to NMM, but the extra powers of Vander¬ 
monde determinants are canceled by the HCIZ integral. 

Indeed, it is not difficult to notice that the H2MM is very similar to NMM if we take 
in the former I\ = K. C = tk ■ k = 1,2,... and compare the eigenvalue representations 
( |1.2| ) and ( |1 . 4| ) . We see that the difference is only in the fact that in ( |1.2| ) the eigenvalues of 
two matrices are complex conjugate, whereas in (|1.3|) they form two independent real sets. 
The number of the integration variables is the same, only the contours of integration are 
different. It suggests that for a large class of potentials W with good convergence properties 
at infinity (in particular, for many polynomial potentials) the partition functions should be 
equal, and the correlators of the quantities Tr <f> n , Tr <f>f m should be the same as correlators 
of TrX n , Tr V'". In what follows we will mostly consider the symmetric case K — K and 
in particular examples even restrict ourselves to the case of real coefficients tk — th¬ 
in the large N limit the models will be also equivalent for most of the potentials W. 
The solution to the saddle point equation for H2MM gives rise in general to two sets of 
complex conjugated eigenvalues: x m — z m , y m — z m m = 1,..., N, since it is the only 
way to make the result for the partition function real for a general set of complex couplings 
tk- Hence they will form the same spots of the two dimensional Coulomb charges with the 
uniform density p(z 7 z ) = —ddV(z, z) = 1 as in the case of NMM |24|], 0. 

The simplest demonstration of this equivalence comes from the direct calculation of 
the Gaussian H2MM integral 


z Gauss = V X VY e 


4n \ 


Tr (- XY+t 1 X+t 1 Y+t 2 X 2 +t 2 Y 2 ) 

2 (+ t\t2 + f 1^2 


N 2 /2 


1 - 4 t 2 t 2 J 


exp N 


1 - 4 t 2 t 2 


(1.5) 


which of course coincides with the partition function of NMM with the same quadratic 
potential W(z). The latter has in the large N approximation the distribution of the 
eigenvalues in the shape of ellipse [24], [25[, and the coincidence of results confirms our 


statement that in the saddle point approximation the eigenvalues Xk,Vk will be complex 
conjugate and will also form the same ellipse. Both NMM [j9[] and H2MM [|10[ were proposed 
as matrix models describing the two dimensional Laplacian growth processes. 

In the next section we will demonstrate the planar diagram technique for the one 
and two cut H2MM, relating it to the Ising model and to the two phase Ising model on 
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planar graphs, respectively. In section 3 we will reproduce (by an unusual method) the 
solution of the two matrix model in the planar approximation. Using this solution we will 
build in section 4 the general algebraic curve (in general not hyperelliptic) describing the 
multi-support two matrix model. We will describe the topology of its Riemann surface 
and its possible degenerations into lower genera. In section 5 the free energy of the model 
will be presented as a tau-function in terms of the variables corresponding to the periods 
of holomorphic differentials of the curve. The cubic case with the two cut degeneration 
will be considered in detail and a rather explicit solution for its free energy will be given. 
In section 6 we will describe the connection of the two matrix models (and of some of their 
generalizations) to the calculation of effective superpotentials of J\f — 1 super Yang-Mills 
theory with two adjoint chiral multiplets and an appropriate tree superpotential, in relation 
to the conjecture of (T9|] . In section 7 we will sketch out the construction of the algebraic 
curves for a very general class of solvable matrix models with the tree-like interactions of 
the matrices. 


2. Combinatorics of planar graphs of the multi-support two-matrix model 


The equivalence of NMM and H2MM is useful to give a combinatorial interpretation 
to the NMM in terms of the planar graph counting. The NMM does not have such a direct 
interpretation but the H2MM does have it. Indeed, the equivalence between ( |1.1|) and ( |1.3| ) 
suggests the following recipe: if we want to calculate the partition function or the correla¬ 
tors of traces of normal matrices (without mixture of two matrices inside each trace!) we 
simply have to calculate the corresponding quantities in the corresponding H2MM, taking 
arbitrary hermitian matrices X. Y instead of the commuting complex matrices ( I>, <f>t, and 
the same set of complex conjugate couplings t n . t n . We can do it by all available methods 
in the H2MM: saddle point approximation, orthogonal polynomials or loop equations, or 
even perturbatively in the couplings, by the direct planar graph expansion. 


2.1. One-support case: Ising spins on planar Feynman graphs 

Let us remind that the H2MM was used in [|3;] to define and solve exactly the Ising 
model on dynamical planar lattices. The role of these lattices is played by planar Feynman 
graphs and the positions of each spin correspond to two types of interaction vertices ( X- 
vertex and Y-vertex as spin-up and spin-down). In this respect, we can say that the NMM 
also describes the Ising model on planar graphs. The phases of the complex couplings 
correspond here to some generalized (imaginary) magnetic fields, whose values depend on 
the phases of couplings. 
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For example the two matrix model with the cubic interactions describes the statistics 
of Ising spins on $ 3 -type planar graphs (or, due to the Kramers-Wannier duality, on planar 
triangulations): 


crH) 

Z'N 


[ 7 ,A,A]= f 


( 2 , 1 ) 


It corresponds to the following choice of couplings in ( |1.1|) : t\ — t\ — 0, t 2 — t 2 — y/2, 
£3 = A/3, £3 = A/3. Note that this choice of couplings does not lead to the loss of generality: 
any cubic potential can be brought to that used in (0) by constant shifts of matrices 
X —» X + const, Y —► Y + const and the phase rotations X —> e l9 X , Y —> e~ %e Y. 

Let us denote A — g e iH . The planar $ 3 -graphs are generated in the large N limit by 
the expansion of ( |2.1| ) in powers of g. and the corresponding Feynman rules can be given 
the statistical mechanical interpretation in terms of the Ising model on ( I )3 -iype planar 
graphs with the temperature an d the imaginary constant magnetic held iH |§. 

The solution of two matrix model describing the Ising spins on planar graphs corre¬ 
sponds to the situation when the eigenvalues for both matrices form one connected support 
around the classical minimum of the potential corresponding to $ = <fd = 0 . 


2.2. Multi-support case: multi-Ising phases 

Both NMM and H2MM admit in the large N limit the multi-support solutions, in 
analogy with the hermitian one matrix model, where they were studied from the point of 
view of their relation to the hyperelliptic curve in the works [|I7j , [ 2 ( 


m. and recently in 


m and j20| • In the case of NMM the eigenvalues Zk, Zk are distributed with the constant 


density in a set of disconnected spots on the complex plane z. Our main purpose in this 
paper is to describe and classify such solutions from the point of view of the underlying 
algebraic curves. 

We will work in the notations corresponding to the complex conjugated Zk , Zk of the 
NMM. However, all results will be true for the independent Zk,Zk, as in the H2MM. 

The eigenvalue supports appear around the extrema of the potentials. Note that in the 
sense of analytic continuation one can also formally “fill up” all extrema of the potential 
and not only the minima. This leads to more general solutions |l£j leading to important 
physical applications. The extrema of the potential V(z, z) — —zz + W(z) + W{z) are at 
the points defined by the system of equations 


z = W'(z), z = W(z ) 


( 2 . 2 ) 


In general, for the potentials of degree (n + 1) we have n 2 extrema. 
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To be more concrete let us study the case of cubic polynomial potential (related to 
the one in ( | 2 . 1 |) by a simple shift of variables): 


V(z, z) = -zz + T(z + z) + | (z 3 + z 3 ) 

o 


(2.3) 


with real couplings T and g, and fill out only the two extrema obeying the reality condition 
z — z. The classical equations 
of a ” classical” curve a 


|) for the extrema of (|2.3|) can be rewritten in the form 

1 


3 ( 


9 


z-gz 2 - T ) (z - gz 2 


T) = 


= z 2 z 2 - - (z 3 + z 3 ) + — (z 2 + z 2 ) + \zz - (z + z) + ^ = 0 

9 g g 2 g 2 g 2 


(2.4) 


and the solution has n 2 — 4 extrema, two extrema for z — z and another two for z + z — 
— 1/g. The potential can be expanded around the extrema z — z as follows 


777 

V(Z,Z) = ~(Z - Z a )(z-Z a ) + —^ ((Z - Za ) 2 + (Z- Z a ) 2 ) + 


where a — 1 , 2 , 


+ | ((z - z a ) 3 + {z- z a ) 3 ) ± const 


1 


h,2 = ^(1 ± V 1 - 4T g) 


(2.5) 


( 2 . 6 ) 


and rai^ = 2 ^ 1 , 2 . We will not consider the filling of the spots corresponding to other 
extrema, with z + z — —1/g (later we will discuss this fact in a more general context). 
Let us regroup the eigenvalues into two groups and denote: 


(z i - z 1: ..., z Nl - z x ) = (ai, ...,a Nl ) 

(zni.+i — z 2 , ■ ■ ■, zn — z 2 ) = (bi ,..., 6at 2 ), Ni + N 2 — N 


(2.7) 


and the same for the conjugated variables!, corresponding to their positions in the first or 
second spot, respectively. Now we can use the eigenvalue representation ( |1 . 2|) and rewrite 
this integral in terms of hermitian matrices A, A , having the size iVi x Aq and B, B having 
the size N 2 x N 2 and a pair of complex rectangular anticommuting ghost matrices C, C 
having the size N\ x N 2 , as follows: 


Z N [t,t } = / VAVAVB VBVC VC e NTrS{A ’ A ' B ^' C ' d) 


( 2 . 8 ) 


4 Compare it to the curve for the one-support solution in [27] 

5 In the case of the H2MM one can imagine the situation when the filling numbers of the 
variables z % are not the same as for z% s: N[ A Ni, N 2 A W; we don’t consider this situation here. 
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Fig. 1: A planar graph of the two matrix model with two eigenvalue supports. 
There are two phases here: thin line phase (inside the circles) and thick line phase 
(outside the circles). Each phase corresponds to two different kinds of Ising spins 
having a different temperature: the spins looking “up” are located in the triple 
vertices made of solid double lines, and the spins looking “down” are located in 
the triple vertices made of dotted double lines. The three types of propagators 
inside each phase (solid, dotted and mixed) describe the interactions depending on 
the mutual orientation of the neighboring spins. Along the interphase (ghost) lines 
drown by circles, the spins have the same orientation. 


where 

S(A, A, B, B, C,C)=-AA+ ^-(A 2 + A 2 ) + f (A 3 + I 3 ) 

z o 

- BB + — (B 2 + B 2 ) + -{B 3 + B 3 ) (2.9) 

2 3 

- mC ] C - m&C - gC ] CA + gCC ] B - gC^CA + gCCA B 


and to = y/1 — ATg. To rewrite the eigenvalue integrals as the matrix ones, we used the 
HCIZ formula for the U(Ni^) group integral 


f Mm*,) cc 


det t? e aiaj 
A(a)A(a) 


and similarly for b , b. The matrices C , C served to exponentiate the cross-products 
life m( a k —bm)(dk — b rn ) in the Vandermonde determinants in ( |1.2|) (see the similar method 
for the one matrix model in 0,0). 

Now we can give the model ( [2.9|) a combinatorial interpretation in terms of the planar 
graph expansion. Namely, we expand Q2.8|) in the cubic coupling g. for the fixed mi, to 2 , to. 
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The diagram technique consists of the following elements: 


propagators : 


/ ~ ~\ 2toi / . ~\ 2 


<^>o = {AA) q = 
(BB) 0 = (BB)_ = 


(A A) = 


0 \ ^ A / 2 i 5 \ -* ■*-* A / 2 i’ 

\ / o mf — 1 \ / o — 1 

2m2 / „ 2 


2 i ? 

0 777,2 — 1 

<e f e>„ = (e f c) = i/m 


BB) = 


0 m 2 — 1 


( 2 . 10 ) 






vertices : (AAA) 0 = (AAAj = ( BBB) 0 — 


BBB 


- (C'CA) 0 = (CC'B) 0 = - {&CA) o = (<?<5»b) o = 


( 2 . 11 ) 


Each type of the propagators (C'^C') Q and (^C^C^j forms closed loops on Feynman graphs, 
each loop entering with the factor (—1). 

A typical planar graph for the two-support model is presented on fig. 1. Let us classify 
the index loops of each planar graph (or a graph of a fixed topology) as carrying the index 
i = 1,, Ni (solid line), or the index %' = 1,..., Ni (dotted line). The ghost loops (drawn 
by a double line formed by a thick and a thin line) will separate two phases on the planar 
graph: one described by the matrices A, A (thick line phase) and another described by 
the matrices B,B (thin line phase). Each of these phases corresponds to the dynamics 
Ising spins of 1-st and 2-nd kind, as described in the previous subsection for the single 
support case. At the phase boundaries formed by the ghost loops, the spins have the same 
orientation (two types of ghosts C and C correspond to two possible orientations). Each 
solid index loop contributes a factor Ad, and each dotted index loop - a factor N 2 . 

Let us note at this point that for gT < 1/4 we always have real mi > 0, m 2 > 0, but 
the determinants of second derivatives of the action at two different extrema are m\ — 1 > 0 
and mj — 1 < 0 correspondingly, which means that the first extremum is the true minimum, 
and the second is a saddle point of the potential. 

The last comment about this diagram technique: as it was done in [|2(J for the two- 
cut one matrix model, we can do the following formal operation with each graph, without 
changing its contribution: we can change the contribution of each ghost loop from (—1) to 
1, change the sign of each BB , BB and BB propagator, and the sign of N 2 I. 

All this means that we can consider instead of the matrix model (|2.8|) -( |2T| ), the matrix 
model with the action: 

S(A, A, B, B, C,C)=-AA+ ^(A 2 + A 2 ) + |(A 3 + A 3 ) 

Z u 


+ BB 


f(B 2 + B 2 ) 


+ | (B 3 + B 3 ) 

o 


- mC^C - m&C + gC^CA + gCC^B + gC ] CA + gCC ] B 

( 2 . 12 ) 


6 we will see in the next sections that gNi/N and gN 2 /N can be viewed as independent 
variables in the planar limit, as in the one matrix model case 








Here C and C are already the usual commuting complex N± x N 2 rectangular matrices 
(the sign of N 2 is again normal here). We also changed the variables as follows A —> — A , 
A -> -A. 


In this new representation of the same model, the perturbative (/-expansion goes 
around the true minima of the potential, and the contributions of planar graphs are pos¬ 
itive. The planar expansion of this matrix model define the statistical mechanical model 
on random dynamical graphs describing a two phase system, each phase corresponding to 
the system of Ising spins with the ferromagnetic boundary condition on the phase bound¬ 
ary. Note that since we have two independent “cosmological constants”: the coupling g 
and N 2 /N. but only one parameter T related to the (different) temperatures of two kinds 
of Ising spins. Hence we cannot make both types of the Ising spins critical at the same 
time. We need for that higher powers of the potential. In a sense, our multi-cut solution 
generalizes the ADE models proposed in |28|] , |2^ . 

Let us conclude this section by noticing that much of what we did here on the two 
support case can be carried over to the 4-support case of this model and to the multiple 
supports for the potentials of higher degree. However, unlike the cubic case with real 
couplings, the details are difficult to work out. Below we will return to the cubic potential 
and discuss in detail the generic 4-support- structure. We will also see that the generic 
4-support solution has a very natural 2-support ” degeneration”, corresponding precisely 
to the perturbation theory considered in this section. 


3. Solution of the model in the planar limit 


Let us now turn to the solution of the two matrix model in the planar or large N 
limit. As is well-known in this case the computation of matrix integrals ([LT|) or ED can 
be reduced to the solution of the saddle point equation. 

The saddle point equation for the model ( |1.2|) with the eigenvalues z\..... zn (or 
analytically continued saddle point equation for (|1.4|) ) reads 

1 


Zk = 


W'(z k ) + ^ 




z k - z j 


(3.1) 


together with the complex conjugated equation. For the resolvents of distributions of the 
eigenvalues 

G(*) = ft / Hr — LA , <3(ft) = ft / Tr (3.2) 


it can be written as: 


z - $/ 

z = W\z) + G(z) 
z = W\z) + G(z ) 




(3.3) 
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where the resolvent has the usual asymptotics at large z or z for the finite supports: 


G(z) —> to/z + 0(1/z 1 ), G(z) —> to/z + 0(1/z 2 ) 


(3.4) 


To hx the resolvents in ( |3.3|) we have to impose the condition that the functions z(z) and 
z(z) are mutually inverse: 

z(z(x)) — x (3-5) 

To justify this condition we recall once again that the solutions of these equations describe 
the spots of Coulomb charges with the uniform distribution of the eigenvalues with co¬ 
ordinates (zi, zi), ..., (. zn , zn) with the density p(z, z) — 1. The boundaries of the spots 
are in general smooth curves in the complex plane z depending on the couplings of the 
potential. To hx the form of these boundaries it is enough to consider the eqs. ( |3.3|) at 
the boundary. Then both equations should define the same curve z(z). It means that the 
solutions of these equations, z(z) and z(z) respectively, should be mutually inverse, i.e. 
obey eq. (|5~5|)B. Note that in general z should be treated as an independent function on 
the complex manifold (with involution) and it becomes literally complex conjugated to the 
function z only on some real section - the real analytic curve in the sense of 0, which 
is just a boundary of the eigenvalue distribution. To avoid further misunderstanding in 
what follows we will denote this function as z(z), so that z(z) — z (i.e. is literally complex 
conjugated only on the boundaries of the spots). 

In the quasiclassical, or dispersionless limit one considers the free energy of the matrix 
ensembles to be defined as a ” planar” limit 


jF(t, S) — lim ( Tl 2 log Z 


(3,6) 


implying N —> oo, % —> 0 with NTi — f 0 being fixed. In (|3]6]) t denote the parameters of 
the potential V(z, z) — —zz + W(z) + W(z) while S are the new variables directly related 
to the ” filling numbers” of various eigenvalue supports. More strictly, by the planar limit 
([lb|j one usually understands the solution to the variational problem 


F oc / V(z, z)p(z, z)d z - d zid z 2 p(zi,zi)p(z 2 ,Z 2 )log\zi-Z 2 \ + 


+ 5>(/ p(z,z)d 2 z-S c 


(3.7) 


7 A more rigorous derivation of these equations from the method of orthogonal polynomials 
can be found in || , || . 
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which is the ”stationary phase” condition for corresponding matrix integrals (1.2), (1.4). 
Note that the normalization of the density at different supports is achieved by the Lan- 
gange multipliers v a - The dispersionless tau-function can be obtained from ( |3.7| ) by the 
substitution of the saddle point solution p = p c of the saddle point equation 

W, = 0, 


or 


Sp(z, z ) 

v a = d 2 z'p(z',z')log\z a -z'\-V(z a ,z a ) 


(3.8) 


for any point P a — (z a , z a ) belonging to one of the supports, labelled by a. 

In the next two sections we will first discuss the structure of the complex curve 


for two matrix model and then define the free energy ( |3.7|) as a (logarithm) of quasiclassical 
tau-function. 


4. Complex curve for the two-matrix model 


The reality condition suggests the following ansatz for the solution to ( |3.3|) : 

F(z, z) = fcj zl z 3 = 0, (4.1) 


i,3 


with the coefficients obeying the symmetry: f, t j — fj t □. Due to this symmetry the equation 
( |T5|) will be automatically satished. 

The coefficients /,; 7 - can be partially hxed by the asymptotics: 


2WUoo= i W , 'W + 7 + o(^ 


following from (|3.3|) and (3.4), but the number of parameters of the potential grows linearly 
with its degree while the number of coefficients of ( [4. 1| ) grows quadratically. The rest of the 
parameters will correspond to the eigenvalue filling numbers for various spots (supports of 
the eigenvalues on the 2 plane). Altogether they will play the role of moduli of complex 
structure of the algebraic curve defined by the eq. ( |4.1| ). 

One may also think of the analytic curve ( |4. 1| ) as of the algebraic form of the large 
N loop equations in H2MM (but not in NMM!) for the resolvent G(z) |]30H , fl27H , or for the 
matrix model in external held [pT|, where it was first proposed. 


4-1. Structure of the curve 

We can precise the algebraic equation of the curve (4T) for the mutually complex 
conjugated potentials of a degree K — n + 1 (with a few explicitly given highest degree 
terms) 

PV,, ^ ^ n~n _j_ „^n +1 J_ 7 ,~n+1 j_ E Jiy 

i,je(N.p.)+ 


F(z, z) = z n z n + az n+1 + az n+1 + V f ij z i z j - 0 (4.2) 


8 since any such curve should be consistent with its real section, the eq. F(x + iy,x — iy) — 
V(x,y) =0 with real coefficients. 
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Fig. 2: The Newton polygon for the curve (4.2). The highest degree terms in 
«B determine the shape of the polygon and the integer dots inside it count the 
number of holomorphic differentials, or genus of the curve. Clearly this number is 
equal to the area of ’’dual” square except for one (black) point, so that g — n 2 — 1. 


where the first three terms correspond to the three points on the boundary lines of the 
Newton polygon (square in this case) on fig. 2, and the sum over ( N.P.) + in the last term 
means the sum over the points inside the Newton polygon (including the points on both 
axis not marked on fig. 2). For example, there are 8 terms in this sum for n — 2. 

One may compare this equation with ( |2.4|) and see that the higher degree terms are 
always fixed by the ” classical” equations on extrema of the matrix model potential. The 
properties of the curve ((O)) can always be easily established via the Newton polygon on 

fig- 2. 

Counting the number of integer points inside the polygon one finds that the number 
of holomorphic differentials, or genus of the curve, is equal to 


9 = n 


1 . 


A simple basis for the holomorphic differentials can be chosen as 


dz 


dv,,j — z l z 3 „ 
13 F. 


—z Z J 


dz 


(4.3) 


(4,4) 


with the degrees i — i' — 1 and j — j' — 1, where ( i',j ') e N.P. are coordinates of the 
points strictly inside the Newton polygon, without the boundary points (see fig. 2) i. 


9 For example, for n = 2 there are three points inside the polygon: i',j' > 0 and i + j' < 2, 
then the holomorphic differentials are labeled by i, j > 0 and i + j < 1. 
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Finally let us point out here that for the models with non-symmetric potentials W 
and W one may write the spectral curve equation in a similar way, but it would not obey 
such symmetric properties. For the potentials W and W of degrees n + 1 and n + 1, 
correspondingly, the highest terms will always be of a particular form (one may equally 
use here x and y instead of z and 5 , to show its direct relation to ( |1.4|) ) 

F(z, z) = z n z h + Az n+1 + Bz™ +1 + (4-5) 

(i,j)e(JV.p.)+ 

and the genus of the curve ( |4.5|) is nn — 1. It means that the (n + 1) x (n + 1) square 
on fig. 2 should be replaced by the rectangle of the size (n + 1) x (h + 1) with all other 
elements of the construction remaining intact. Of course, for n — 1 integrating over the 
matrix with Gaussian potential one returns to the 1MM with the (hyperelliptic) curve of 
genus n — 1 (see section 5.1 below). 


4-2. The cubic example 

To understand better the structure of the curve let us first discuss in detail the cubic 
example. Writing eq. ([1^) first with arbitrary coefficients 

F(z, z) = z 2 z 2 + az 3 + az 3 + bz 2 z + bzz 2 + cz 2 + cz 2 + fzz + qz + qz + h = 0 (4.6) 


one finds that its structure is in fact very similar to the equation of the ” classical” curve 
( EOj) . Indeed, the eq. ( f4.6|) should be consistent with the asymptotics 


= W'{z) + G(z) = kt kZ k ~ l + O (^) = A^ 2 + 7 z + 77 + O (- 


k= 1 


(4.7) 


Substituting ( [4.7| ) into the eq. ( |4.fi| ) and collecting the coefficients in front of the terms z 6 , 
z 5 and z A one gets 


1 _ i 
A g 

b =h° 

c — IL _ 21 2 — = — 

A AA A 2 g 


(4.8) 


and their complex conjugated counterparts, i.e. the coefficients at higher degree terms are 
indeed completely fixed by parameters of the potential ( |2.3|) . Four lower degree coefficients 
/, q. q and h correspond to the bipole differential and three holomorphic differentials 0. 
Their classical ” expectation values” are presented in eq. 0). 


10 Note, that the curve ( |4.2D , ( |4.6| ) is written implying some reality condition onto the coeffi¬ 
cients, but as usual, the deformations of these coefficients should be considered as independent 
complex variables. 
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Fig. 3: Cubic curve as a cover of z-plane. 


Let us now present the curve ( |4.(j| ) as a Riemann surface of a multi-valued function 
2 ( 2 ). Then it can be thought as a three-sheet cover of the complex 2 -plane. On the first, 
physical sheet, there are no branch cuts at 2 —> 00 , as it follows from the asymptotics ( |4.7| ) . 
This asymptotics should be supplemented by the ”complex-conjugated” asymptotics 


2 = W'{z)+0[ l 

Z 


(4.9) 


on ’’unphysical” sheets. Then it is clear from (T7) and (O) that on the physical sheet at 
infinity 2 cx 2 2 , while on two unphysical sheets 2 cx yfz and two infinities on unphysical 
sheets are ’’glued” by a cut. 

The branch points at 2 -plane are determined by zeroes of the differential dz, or by 
F~ — 0. Considering the simplest non-degenerate case of the curve ( |4.0| ) 


2 2 2 2 + az 6 -t- az 6 + h — 0 


(4.10) 


it is easy to see that there are nine branch points in the 2 -plane without infinity z — 00 
(of course, one comes to the same conclusion looking at the Cardano formula, or from the 
index theorem, see below). 

The structure of the curve can be then presented as on fig. 3. It is clear from this 
picture that the curve can be presented as two copies of P 1 ’’glued” by four cuts, i.e. 
in general position it has the genus g — 3. There are two ’’infinities” 2 = 00 , 2 = 00 , 
one of them is a branch point. We have shown schematically the possible cuts and the 
corresponding choice of canonical A-cycles. In the classical situation (|2.4| ) one has two 
parabolas intersecting at four points, and under quantum resolution these points turn into 
four cuts connecting two P^s in fig. 3. 
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Fig. 4: Generic curve of the two matrix model as a cover of the z-plane. In contrast 
to fig. 3 each fat line consists of a stack of n cuts. 


To conclude the picture, let us make few comments about the curve ( |4.2| ) for a potential 
of a generic degree n, i.e. when W'(z) ~ z n + .... This curve (see fig. 4) can be again 
presented as two P lJ s glued by n stacks of cuts. One of these P 1 ’s corresponds to the 
’’physical sheet”, the other one is glued at the oo_ from n copies of ’’unphysical” ^-sheets. 
Each stack consists of n cuts, so their total number is n 2 among which one can choose 
n 2 — 1 independent, whose number is equal to the genus of this Riemann surface. 

The differential dz has always a pole of the second order at oo + on the upper, or 
’’physical” sheet, and a pole of the order n + 1 at oo_ since z\ 00 _ oc z n + ■ ■ ■■ It gives 
altogether n + 3 poles and from the Riemann-Roch theorem one concludes that the number 
of branching points, or zeroes of dz is equal to 

#(dz = 0) = n + 3 + 2(n 2 - 1) - 2 = 2n 2 + n - 1 (4.11) 

reproducing nine for n — 2. In general position this gives exactly 2 n 2 branch points, 
producing simple cuts and (n — 1) ramification points, connected by cuts with oo_. 

Finally, let us point out that the structure of the curve (|4.2|) and fig. 4 is consistent 
with the structure of the ’’double” in z and z variables, explicitly seen in the one-support, 
solution H and proposed to be the feature of the multi-support solutions by Krichever 
Indeed, the eq. ( |4.2|) is ’’symmetric” with respect to z and 5 variables, so instead 
of the picture on fig. 4, one can draw a ’’dual” picture of a (n + 1)-sheet cover of the 
5-plane. These dual pictures can be combined together as on fig. 5. Cutting ’’physical” 
sheets from both pictures one may glue them together, forming a double with involution 
z <-► 5. The only delicate point is that these z and 5 sheets should be glued together along 
the boundaries of the spots where z(z) — z and vice versa, in contrast to the picture of 
fig. 4, where the sheets are glued along the cuts on Riemann surface of the multi valued 
function z(z) defined by solution to eq. ([T2|). Both sheets of the double (the lower picture 
at fig. 5) generally have n 2 spots. 
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z 


Fig. 5: Generic curve of the two matrix model as a double of z and z planes. 
One takes the Riemann surface of the function z(z), as on fig. 4, and its ’’mirror” 
Riemann surface of the function z(z) which possesses the same structure. Cutting 
the physical sheets one may glue them together along the (real) curves z — z(z). 



Fig. 6: The boundary of the spot 7 and a cut of a multi-valued function z(z) inside 
the spot. On 7 one has an equality z — z(z) but this is not true on the cut. 


The difference between the boundary of a spot and a cut of a multi-valued function is 
demonstrated on fig. 6. One obviously gets the following relations for the two-dimensional 
and contour integrals 



16 







Fig. 7: The curve (T6) as double cover of the torus. When the torus ( 4.1^ ) de¬ 
generates, the genus g — n 2 — 1 = 3 curve ( ]f(\ ) degenerates into the curve of 

^7red — U 1 — 1* 


which clarify the equivalence of two pictures on fig. 5. The relations ( |4. 1 2| ) allow to endow 
the complex curve ( [4.2| ) (or ( |4.5| ) in the asymmetric case) with a meromorphic generating 
differential zdz. 


4-3. Degenerate curves 


Let us now discuss how the curve ( f4.2|) can be degenerated. The (smooth) genus 


g — n 2 — 1 of the curve ( |4.2| ) decreases if there exists nontrivial solution to the system of 
equations 

F(z,z) — 0, dF — 0 (4-13) 


It imposes certain constraints to the coefficients of fij of the equation (|4.2|), which can 


be found, say computing the resultant of the equations ( ^.13|) . However, these constraints 
cannot be really resolved in a general position. 

To get an idea how the curve (fO|) can be degenerated consider first the cubic case 
( fOl) and let us put all coefficients of this equation to be real. Then, it is easy to see that 
it can be rewritten in the form 


1 


Y 2 + aX ‘‘ + cX 2 + qX + h - X ((3a - b)X + 2c - }f = Y 2 + P(X) = 0 


(4.14) 


where 


X — z + z, Y — zz 


4 ((3« 


b)X + 2 c-f) 


(4.15) 


One may ’’tune” for simplicity the coefficients of the potential ( |4.8|) to get 3a — b and 

can be presented as a double cover 


2 c — f. The formulas ( 4.15 ) show that our curve 
of the torus (|4.14|) with four branch points which are solutions to the equation ( |4.0| ) under 
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Fig. 8: The Newton polygon for the curve ( 4.17 ) gives the genus g * = 


n(n— 1) 


the substitution z — z, where the transformation ( |4.15Q becomes singular. Hence, the 
curve Q4.6|) can be also presented (in addition to fig. 3) as two tori glued by two cuts (see 
fig- 7). 

Now it becomes clear how this picture can be degenerated. Rewriting equations ( [4.13Q 
as 

F s =zF y + P'(X) 

4.16 

F~ z = zF y + P\X) 

one immediately finds that they lead either to z — z or to Fy — 0 and P'(X ) = 0. In the 
second case the torus ( [4.14|) degenerates, while 2 = 5 leads to degeneration of the cover of 
this torus. We will be more interested in the degeneration of the torus since, for example, 
it corresponds to filling of the ” correct” vacua (real eigenvalues) in the perturbative picture 
considered in section 2. 


When the torus degenerates into a rational curve, one gets the Riemann surface ( fl.6| ) 
presented as a double cover of this rational curve with two cuts, i.e. as a Riemann surface 
of genus g — 1 with smooth handles of tori degenerated into (a pair of) singular points 
(see fig. 7). 

The equations of degeneration of the torus ([4 .1 4|) can be easily written using the 
conditions for the double root of the polynomial P(X). Explicitly these conditions acquire 
the form of the discriminant of P(X) or the resultant of the two polynomials P(X ) and 
P\X). 

Now, in the general case ( |4.2|) with real coefficients the substitution analogous to 
( [4 .1 5|) brings it to the form 

Y n + X n+ 1 + ... = 0 (4.17) 
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Fig. 9: The general curve ( U-$) as a double cover of the curve ( 4.17 ) with a genus 
g* — n(Ti 2 ~ 1) . Similarly to fig. 7, when the curve ( 4.111 ) completely degenerates into 
a rational curve the curve of two matrix model ( \4.2 | ) degenerates into the curve of 
genus g red = n - 1. 


where by dots we denoted monomials of lower degrees in X and Y , and there are no ” mixed” 
terms in the eq. ( |4.17| ). The genus of the curve ( |4. 1 7| ) can be again easily computed by 
the Newton polygon, which gives 


9 * = 


n(n — 1) 
2 


(4.18) 


In the same way one may present the generic curve of the two matrix model ( |4.2| ) as a 


double cover of ( |4 .17|) with 2 n branch points. Indeed, the Riemann-Hurwitz formula 


2-2g = #S-(2- 2<7q) - # B.P. 


(4.19) 


where # S is number of sheets and B.P. is number of branch points, gives for g — n 2 — 1 
and go — g* exactly B.P. — 2n. It means that the generic curve of the two matrix 


model Q4.2|) can be presented as a double cover of the curve ( |4.17| ) with n cuts, and when 


the curve ( 4.17 ) degenerates into a rational one, the curve (T2) has the genus 


fi^ed — Tl 1 


(4.20) 
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4-4- Rational degenerations 


Let us finally say a few words about the rational degenerations of (L2), i.e. when its 
(smooth) genus vanishes. A particular example of such a totally degenerate curve is given 
by the ”classical” curve ( |2.4| ). but the rational case can be easily studied for the generic 
values of coefficients in (|4.2|) , i.e. without any reality restriction. 

In such situation eq. ( |4.2|) can be resolved via the (generalized) conformal map 


— rw + 


fc=o 


Uk 

w k 


(4.21) 


/ - k 

= — + > u k w 
w 

k=0 


and the substitution of ( |4. 2 lp into ( |4.2| ) gives a system of equations, expressing all coeffi¬ 
cients fij in terms of parameters of the conformal map ( [1.211) . 

Indeed, substituting ( [4.21|) into (|4.2|) and computing the residues one finds that the 
expressions 

i?z[F]=res^— w l F(z(w),z(w))^j =0 (4.22) 

for l — —n(n + 1 ),..., n(n + 1) form a triangular system of equations onto the coefficients 
fij. It means that each of the equations ( |4.22| ) is linear in one of the coefficients, and can 
be resolved step by step, starting from the ends of the chain. 

For the cubic potential [n — 2) the solution is 


a — 


r 

u 2 


b = 


U\ r 

u 2 


2 uq 


u i ruo 
U2 


/ = r- 2« 2 m 2 +4m 0 mo + 


+ Uq — 2 r ui + 3 


r 2 uq 
u 2 


r 3 ui 
U2 U2 


r 


q — —3 


r u, 


o 


U2 U2 
-2 


tuiUq r tti ut rw 0 ui 

— Ul U\ — 2 —3 - 1 - 3 - 2 - 


+ 2U2U2UQ — U2U 1 


U2 

+4 Uo r Ul — 2 Uq Uq — 


r 2 ul 
U2 


U2 

4 - 
r uo 

U2 U2 


U2 u 2 U2 

3rU,2Ui + Uq Ui Ui 


2 ruiui 2 2 - , n r u i u o u o , „ r 3 «i 


u 2 


- r u 0 + 2 


u 2 


+ 3 


u 2 


TUI Uq 
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u 2 u 2 


u 2 u 2 
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—u 1 Ui U2 U2 — 3 r U2 U2 + 3 U 2 T Ui Uo — 2 U2 U2 Uq Uq ~ U± r Ui + 3 Uo r U2 U\ 


together with the ”complex conjugated” expressions for a, b. c and q. where one should 
replace Uk by Uk and vice versa. Resolving (|4.22| ) one gets the explicit description of the 
rational degeneration of the curve ( |4.2| ) in terms of the coefficients of conformal map ( |4.2 1 | ) . 


However, in general situation they are only implicitly dehned through the parameters of 
the potential V(z, z). 


5. Quasiclassical tau-function 


Let us now define the partition function for the two matrix model ( |3.7|) in terms of 
the quasiclassical tau-function introduced in |T6[]. First, we discuss a simpler example of 
the one matrix model and then turn to the particular features of the two-nratrix case. The 
hyperelliptic curve of the one matrix model was first discussed in [17|, and recently, in the 
most general form including all the extrema of the potential, in |19[ (see also ||20|| and ||33||). 


5.1. One-matrix model 

The complex curve of the one matrix model 


Z= d<te T ' w "W 


(5.1) 


with the model potential 


«*) = X >*** _1 (5.2) 

k =1 

comes from the very simple loop equation G 2 + 2W^(X)G — /(A) = 0 (see, for example, 


131 


or [|J). It is always hyperelliptic, i.e. can be rewritten in the form 

y 2 = VF'(A ) 2 + /(A) 


(5.3) 
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with y — G + W' n and the moduli hidden in the coefficients of the polynomial 


71—1 

/(A) = X>A*. (5,4) 

k=0 

The generating differential is chosen as 

dS 1MM = — ydX (5.5) 

2 t n y ’ 

and additional variables, corresponding to the eigenvalue filling numbers, can be introduced 
through its periods 

S t = j dS 1MM (5.6) 

directly related to the integrals of density over the eigenvalue supports. Then 


OdS 1MM 

dSi 


duji 



(5.7) 


where the derivatives are taken at fixed coefficients {ti} of the potential W^( A) ( |5.2|) . As 
usual, the periods dual to ( |T . fj| ) are given by the integrals over dual cycles 


jjlMM 


dS 


1 MM 


Bi 


(5.8) 


To complete the set of parameters of the model, we have to add to the filling numbers 
( lint) and coefficients of the potential ( |5.2| ) the variable 

reSoo + (dS 1MM ) = -res OQ _ ( dS 1MM ) = = t 0 , (5.9) 

zt n 

so that 


dds 1MM 

dt 0 

The dependence of {fk} with k — 



X n ~ l dX 1 (V aft X l dX 

~ tn ^ + ~ 2 h w ~' 

0 ,1 ,..., n — 2 upon to is fixed by 


71—1 


dX 


1 n—2 


df k X k dX 


2 ^ dt 0 y 

k =o u y 


= 0 


(5.10) 


(5.11) 


which gives for i 
differential ( |5.10|) 


= 1 ,..., n — 1 exactly n — 1 relations on / 0 , fi ,..., f n - 2 - 
can be also rewritten as 


The bipole 


dfl± 


dds 1MM 

dt 0 


/I hit 


f g(P,oo + ) A 
\E(P, oo-)J 


(5.12) 
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where E(P,P') is the Prime form on ( |5.3|) . Obviously, the differentials ( |5.12| ) obey the 
properties 

reSoo + (iO-|- = — res 00+ dCi± — 1 

/ (5-13) 

j) dfl± = 0 , i = 1 ,..., n — 1 

To complete the setup one should also add to ( |5.8|) the following formula 

/» 00 _|_ 

dS 1MM (5.14) 


n n = 


which can be regularized in the usual way presenting the puncture at infinity as a degen¬ 
erate handle. The partition function of the multi-cut solution of the 1MM is defined now 
in terms of the quasiclassical tau-function p 1MM obeying the equations [|T7| , [|19 


= n; 


qjzIMM 

dS} MM 

qjtIMAI 

dtn 


1 MM 


(5.15) 


= n r 


In the papers ||19|| , instead of to the parameter S — t 0 — = S n was used. This is a 


non-standard definition of homology basis on (|5.3| ) and it gives rise to the divergences at 
infinities. However, the basis of |Hj is related to the canonical one by a linear change of 
variables, where no divergences appear (except of the trivial one in ( |5.14|) ) and the inte- 
grability of (|5.15| ) (the symmetry of second derivatives) follow from the Riemann bilinear 
relations, including the symmetry of period matrix of 


5.2. Two-matrix model (general potential) 

In the same way, the filling numbers can be defined for the two matrix model 




2iri 


z —th spot 


dz A dz — - ® zdz = ® dS 

2 ™ J Ai Ja, 


2 MM 


i.e. as periods of the generating differential 

dS 2MM = 


27 li 


zdz 


(5.16) 


(5.17) 


under the appropriate choice {Ai} for the basis of T-cycles on the Riemann surface ( 14.20 , 
(or (|4.5|) ). This is illustrated by fig. 4 (or by its particular cubic case fig. 3), taking into 
account fig. 6 and eqn. ( 4.12|) . From (|5.16|) one still gets in the same way the analogs of 


the formulas ( |5.7|) 


dds 2MM 


dS t 


= du>i 


dujj = Sij 
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(5.18) 
































where the canonical holomorphic differentials (now on the curve (4.2)) are certain linear 
combinations (with moduli dependent coefficients) of g — n 2 — 1 ” lower degree” holomor¬ 
phic differentials ( |4.4D . 

The derivatives of ( |5.17|) with respect to the coefficients of the equation ( |4.2|) can be 
computed in the standard way. Choosing z as a covariantly constant function one writes 
for ( fTj ) 

Fzdz + SF = 0 (5.19) 


where SF = Y Sf, : jZ l z J is a variation of only the coefficients of (|4.2|) . Then the variation 


of ( |5.171) gives rise to 


Szdz = SF^F- = - ^ SfijZ l z J ^ 


(5.20) 


Expression ( |5.20| ) contains a decomposition of the variation of the meromorphic differential 
( [5.17 ) over some basis of meromorphic and holomorphic differentials on the curve (4.2). 
It is easy to check that the coefficients fij corresponding to the meromorphic Abelian 
differentials of the second kind can be expressed through the parameters of the potential 
V(z, z) of the two-matrix model, namely, through the coefficients of its harmonic part 
W(z) + W(z). The corresponding relations follow from the fact that the complex curve 
( |4.2|) should satisfy the asymptotic expansion of the branch 


/l \ n+1 

5 = W'(z) + O - = ^ kt k z k ~ x + 


fc=l 


Oil 

Z 


(5.21) 


which gives rise, e.g. to a~ l — —(n + l)t n +i e l c - 

The rest of coefficients consists of the coefficient corresponding to the third kind 
Abelian differential (with the first-order pole) and the holomorphic differentials Q4.4D . Fix¬ 
ing the coefficients expressed through the parameters of the potential in (|5.20|) and taking 


appropriate linear combinations, one arrives at ( |5.18|) . 


As it will be shown below, the dependence of the free energy Q3.71) upon the filling 


numbers (|5.16|) is defined by 


OF 

dSl 


dS 


2 MM 


Bi 


2tci 


zdz 


(5.22) 


Bi 


where {Bi] are the canonical dual cycles A, o Bj — Sij. The integrability of (|5.22|) follows 


from the symmetry of the period matrix of the curve ( |4.2| ) . 
Now, as in the one matrix case, one should also introduce 


res 


00 _|_ 


(dS 2MM ) = —reSoo_ ( dS 2MM ) = to 


(5.23) 
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adding the bipole differential 


dsl±= 8dS^ = diog f E ( p ,c° + ) 


Otr 


E(P, OO—) 


_ra—1 ~n — 1 
= Z Z 


dz 


dfij , 

holomorphic 


CX/6 \ -V 


(5.24) 


where coefficients ( -^) L are fixed, as in ( |5.11|) , by j> A dfl± = 0. The variables ( |5.16Q 


dto ’ ''E—- v’> J JAi 

should be directly identified with those introduced in (3.7) for a — 1,. .., n 2 — 1 while 

S n 2 = to - El 


S',-. The second kind Abelian differentials are defined as derivatives 


dfifc = 


DdS 2MM 

dt k 


dOjt = 


dds 2MM 

dtk 


over the parameters of the matrix model potentials, and the dependence of the coefficients 
fij upon the coefficients of potential is fixed by j> A dQ},- — f A dQk — 0. 

The formulas ( |5.16| ), ( [T22| ) and (|5.23|) , together with the (regularized) equation 


dP 

dt 0 


r°°+ 
^-= dS 


2 MM 


1 

27 ii 


'*00 + 


zdz 


(5.25) 


define the quasiclassical tau-function [It 


It is clear that this definition coincides in fact with the definition of free energy of 
NMM or H2MM (|37F|). Indeed, using the formula ( |3.7| ) one can easily check that 


dT 

dt k 



1 

27 ii 


res 


(. z k zdz ) 


(5.26) 


and 


dT 

dS~ a 


= v a 


(5.27) 


The quantities ( |5.27| ) are in fact nothing but linear combinations of ( p.22|) and ( |5.25| ) . like 
in the case of the 1MM. It can be derived by carefully treating the logarithmic integral 
f dz A dz log | z — z a |, with z a belonging to one of the eigenvalue supports. Indeed from 
(P75|) one can get for v a — V(z a , z a ) 


J dz A dz log | z — z a 


£ 

spots ’ 


dz A dz\og(z — z a ) + c.c. 


= £ 

boundaries 


log (z — z a )zdz + c.c. 


£ 

boundaries 


log(z — z a )zdz + c.c. 


(5.28) 
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Fig. 10: The integral in $5.28 ) can be transformed to the (linear combinations 
of) the integrals over the cuts of logarithms, which turn into the B-periods of 


dS 2 


= Tr^zdz. 
2m 


The integral in the r.h.s. of ( |5.28| ) is an integral of the multi-valued differential defined on 
a Riemann surface with two sorts of cuts: the cuts of the function z(z) and the cuts of the 
logarithms. The integral over the spots can be transformed into boundary integrals, where 
the boundaries should now include the integrals along the branches of the logarithmic 
cuts (see fig. 10), the corresponding contributions to the two-dimensional integrals over 
spots vanish. The latter ones are combined into the integrals around all spots which can 
be deformed to the infinity, modulo the integrals along the logarithmic cuts. Altogether 
these contributions gives rise to the canonical /i-periods of dS 2MM — zdz. Finally, the 
situation appears to be quite similar to 1MM since applying the procedure of ”transfer” 
of an eigenvalue from one support to another [ [IT)[ | one should leave intact the boundary of 
the spot (since adding an eigenvalue directly to the boundary not only changes the filling 
numbers, but also the shape of the domain, related to the parameters of the potential [p5|). 
Instead, as in one matrix case (where one may neglect this problem due to vanishing of 
the density of eigenvalues at the ” boundaries” of a cut) one should put the points z a to 
the branching points of the complex curve where dz — 0. After that, using (|5.28|) and 


dV(z, z) 


dz 


oc 


dz =0 


cut 


z(z')dz' 
z — z' 


one gets the formula ( |5.22| ) . 

Indeed, let us calculate, for example, the difference by procedure of the 

“eigenvalue transfer” from an endpoint zp — z" of the /3-th cut to an endpoint z a — zJ 
of the Q-t h cut in z-plane. It will be given by difference of the corresponding eigenvalue 
effective actions (see ( |1.2|) , ( |1.4|) ) at the saddle point: 

8T dT = + z n~ z n + _ W ( z »} + _ w(z") + 


dSi dS i 


N 

V - ''' , (z Zm ) ( % Zm ) 

+ ^ og ( Z " - Zm )(~z" - z m y 


(5.29) 
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where z' — z{z '), z" = z(z"), with the function z(z) dehned by the algebraic curve ( fT^ ) 
(or ( |4.5| )) of the two matrix model. 

Passing to the continuum limit and introducing the resolvents G(z) and G(z). as in 
|), we rewrite the last term in (|5.29|) as follows 


c ^ G <"> log 


w - ~_) 

(z" - z) 


+ %^l 6(i) log 


Z' -z 


5) 


where the contours C and C encircle all the eigenvalue supports - the cuts in the z and 
5 planes, respectively (these cuts or their stacks are depicted in fig. 5). It is important 
that, according to the definition of the sums in the last term of ( |5.29|) , the contours do not 
encircle the logarithmic cuts along (z 1 , z") and (z', z") intervals. 

Blowing up the contours C, C we will encircle only the logarithmic cuts (note that 
there are no poles at infinity). Calculating discontinuity along these cuts we reduce the 
contour integrals to the ordinary ones: 

~r 

dz z{z) + dz z(z), (5.30) 

Jz" 


dT dT 


dS„ 0S< 


/ — / i // -// , 

= — 2 2 + 2 Z + 


the potentials being absorbed into the functions z(z) — W\z) + G(z) and z(z) — W'(z) + 
G(z), according to the saddle point equations of the two matrix model. 

Integrating by parts in the last term of ( |5.30|) (note that the last integral after the 
change of variables goes along the unphysical sheets of the curve fig. 4) we finally get the 
integral over the dual B, 3 cycle 


dT dT _ I ~ zdz 

dSi dSj J B .. 2ixi 


(5.31) 


which is equivalent to the equation ( |5.22| ) defining the geometry of the planar limit in two 
matrix model. 


5.3. Explicit form of the two-support solution for the real cubic potential 

In the rest of this section we will study in detail the degenerate case of a cubic potential 
with real coefficients, having only two eigenvalue supports (on the real axis), instead of 
four. The degeneration greatly simplifies the calculations. The period integrals can be 
even rewritten in terms of elliptic integrals. 

Indeed, the degenerate torus ( |4.14f) can be presented as 


Y 2 + P(X) = Y 2 


-(X-X 1 ) 2 (X-X 2 )^0 

9 


(5.32) 
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and the new parameters X\ and X 2 are defined by the discriminant equations 
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g 4<r 
IT 
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3 T 3/ 


nx) = —x^+---,x^+ q + — -^-x + h--—-f =0 


9 


9 


1 / 2T 


9 


m)=-W?-£)* + ( 9+ ^-¥')= o 


r 


.9 


(5.33) 

where we used ( |4.14| ) and ( [4.8| ). This system of eqs. fixes, for example, h which is now not 
independent and can be expressed through q and /. 

Note that we have the expansions 


where 


X 1 = ~ + 5X 1 = -- + 0(6f,5q) 

9 9 

X 2 = -^ + T + SX 2 = -^ + T + 0(Sf,Sq), 
49 4 g 


sf=f- 4 

g 2 

s T 

°q = q + -o 

9 2 


(5.34) 


(5.35) 


are deviations of moduli from their classical values in eq. ( |2.4|) . The “classical” value of 
X\ in ( |5.34j) precisely corresponds to the classical solution of ( |2.4|) , of course different from 
diagonal the z = z. 

Substituting z — z into eq. (|4.G|) one gets the fourth-order equation 


z 2 + -z~- + {- ) + P(2z) 
9 9 2 ' 


z 2 H —z -1— — ) H— (2 z — X{) (2 z — X 2 ) — 0 


9 


9 


9 


(5.36) 

corresponding to the four branch points of two remaining cuts. These cuts can be inter¬ 
preted as ” splittings” of the double zeroes of the ” classical limit” of this equation given 
by (|2.6|). Let us denote the splitting as Zi j2 —■► zf 2 , so that zf and z 2 are four solutions 
to ( |5.36|) . Then the only non-degenerate A- and B- periods on the torus are given by the 
integrals of the differential dS 2MM — ^-zdz between these points 


S = 


pj-p2MM i 

-= n = — 

dS 2ni 


zdz = 


1 

27 li 

1 


zdz 


zdz 


(5.37) 
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where z is related to z via 
3 


T f 

z~z+-(z + z)--+ J - 
2 9 9 2 


H— (z + z — ^i) (z + z — X 2 ) — 0. 
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These relations for the only nontrivial period S ( [">.371 ) 011 this degenerate g re d = 1 curve 
should be supplemented by the relations ( |5.23| ) and (|5.25| ), defining the dependence on 
to- As usual, one may choose instead their linear combinations Si — S and S 2 — to — S, 
corresponding to the filling numbers of the two cuts. 

In this way we formulated the explicit solution of the two support two matrix model 
with the real cubic potential. We point here out again that the integrals ( |5.37|) can be in 
principle calculated in terms of elliptic functions. 


6. Connection with SUSY gauge theories 


Recently it was proposed in |TR] to build the geometries of underlying string theories 
for certain Af = 1 SUSY gauge theories by effectively reducing them to the complex curves 
( |5.3|) of 1MM. The curves of the 1MM belong to the same class as the Seiberg-Witten curves 
of the SU(n ) Af = 2 SUSY gauge theories with n — 1 fundamental matter multiplets. The 
understanding of this proposal directly from the field theory was considerably advanced in 
S, and later in [|36|]. Though the parallels between matrix models and four-dimensional 
supersymmetric gauge theories based on the similarity of their integrable structures were 
noticed much earlier (|13[ , the recent observation of Jl9j] contains a direct conjecture relating 
the super potentials in Af — 1 four dimensional theories to the partition functions of the 
multi-cut solutions, like ones considered in our paper. 


According to the proposal of [19] the effective potentials of gaugino condensates 


Si — / Tr Wa^W^yj in a large class of four-dimensional Af — 1 gauge theories can be 
calculated in terms of the planar limit of matrix integrals. For the Af — 1 theory with 
one adjoint matter multiplet (broken Af = 2 supersymmetric theory) the calculation re¬ 
duces to the large N solution of the one-matrix integral ( |5.1| ), in general having multiple 
cuts, as described in the previous section. When the U(N ) gauge group is broken to 
U(Ni) x U(N 2 ) x ... x U(N k ), with the classical VEV’s of different subgroups located 
at the different extrema of the tree superpotential, the matrix model predicts the values 


of k gaugino condensates Si — ( Tr Wa ' > Wfy ) corresponding to the the vector multiplets 

Wa ^ of the gauge subgroups. The effective potential W e ff(Si,r ) as a function of these 
condensates and the complexified gauge coupling r can be related to the free energy of the 
multi-cut solution as follows 


W eff (S, t) =YNj \ 2nirSi - 


8T{S u ...,S k ) 

dSi 


( 6 . 1 ) 


the logarithmic term being hidden in the second term representing the derivative of the 
matrix model free energy, according to the proposal of [|T9[. The variables Si,, S k 


appear, strictly speaking, only in the planar limit of this matrix model and correspond 


29 
































to the eigenvalue filling numbers S, — hNi oc Ni/N of various classical extrema of of the 
matrix action, giving rise to the dependence of multi-support solutions upon the variables 
jD, ( |5.16|) discussed in our paper. 


For the one-matrix model ( |5 . 1| ) the situation looks to be relatively simple since the 
key observation comes from the fact, that from the coincidence of the matrix model and 
Seiberg-Witten curves it trivially follows that 


q 2 jrlM M 

dSidSj dciidcij 


( 6 . 2 ) 


due to coincidence (after fixing the homology basis) of the period matrices. The Seiberg- 
Witten prepotential [J2j as a function of a different set of variables 


di — 


= 0 A 


dW' n Wf df 

y y 2 / 


is defined as 


dT sw 

dai 


A 


Bi 


dwi w; df 


y 


and it is also a quasiclassical tau-function [11 


y 2 /, 

Note that eq. 


(6.3) 


(6.4) 


0 states only that 


the second derivatives of different functions in different variables coincide, but these 
functions themselves are certainly different quasiclassical tau-functions. Such identifica¬ 
tion became possible first of all since the number of multicut variables for the 1MM is 
n— 1 = rank SU(n), what is precisely the dimension of the moduli space of (the Coulomb 
phase of) the SU (n) gauge theory. 

For the softly broken AT — 4 theory, according to the proposal of [|19|, iF(S i,..., Sk) — 
log Z should be calculated as the planar limit of the matrix integral 


Z — / e 




with $i, $ 2 , considered as simple hermitean N x N matrices. Of course not all matrix 
integrals of this kind are calculable. An important case corresponding to the J\f = 1* 
perturbation C(<f>i, $ 2 , $ 3 ) = m Y ^,=considered in [JlfJ and based on the planar 


solution of this matrix model given in [37 


If only two out of three masses are nonzero and equal, the theory possesses J\f — 2 
supersymmetry and its non-perturbative solution is formulated in terms of the elliptic 
Calogero-Moser system . |39|] • |J3| • whose spectral curve |4C]] covers n times the elliptic 
curve. There seems yet to be no naive and direct relations between the corresponding 
Seiberg-Witten theory and quasiclassical tau-function of two matrix model, considered in 
this paper. However, the structure of its complex curve suggests that certain geometric 
parallels between these two theories are quite possible. 
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Let us discuss a possible place of the two matrix models in the context of the proposal 
of [|T9|]. An obvious and interesting generalization of the J\f — 1 SYM theory with one 
adjoint chiral multiplet is the the case of a few multiplets with direct interactions among 
the fields. If one takes, in the case of two adjoint chiral multiplets X and Y, the tree 
superpotential W tre e — — Tr XY + TrW(X) + Tr W(Y) then the function R(Si, S 2 , ■ ■ ■) 
in the general formula ( |0.1| ) should be chosen as the planar free energy of the two matrix 
model with (in general) multiple supports, being calculated according to the equations 
CT) and (ET22D- 

In this way we establish the relation between the algebraic curve of the two matrix 
model and geometry of the supersymmetric theory with the described class of tree poten¬ 
tials. 

The two matrix models obeys a rich phase structure in the space of its couplings. Its 
critical points correspond to collisions of various singularities on the algebraic curve, like 
the collapse of the B-period of the Seiberg-Witten curve corresponding to the appearence of 
massless monopoles. These critical points were used for the complete classification (within 
the H2MM) of the models of (p, q ) rational matter fields interacting with 2D gravity (see 
for example ||). It would be interesting to study the consequences of this well established 
picture for the phase structure of the underlying J\f — 1 SYM theory with two adjoint 
chiral multiplets. 

Let us also note here that an interesting generalisation of the standard H2MM con¬ 
sidered in this paper is the model describing the perturbed quantum mechanics of the 
inverted oscillator, proposed in |LI], The model is related to the dynamics of windings in 
the compactified 2D string theory. For rational R it can be also described by an algebraic 
curve. 

In contrast to the one matrix model, in the two matrix model the number of multicut 
parameters ( |4.3|) grows exactly as the dimension of the SU (n) group n 2 — 1 = dim SU (n). 
Naively, this would correspond to the total breaking of the SU(n ) gauge group in the 
Seiberg-Witten like context, including even the breaking of the corresponding global sym¬ 
metry, or to a theory with a more complicated gauge/matter structure. 

However, an important particular case of the multicut solution of the two-matrix 
model corresponding to the case of real couplings in the potential and only n — 1 real 
eigenvalue supports, giving the curve with genus (|4.20|) . In this case we have the number 
of extra parameters exactly equal to the genus f/ re <i = n — 1 = rank SU(n). which might 
be more appropriate for the study of the symmetry breaking in the corresponding J\f — 1 
SYM theory. 

There exist a few obvious generalizations of the multi-matrix interactions, giving rise 
to a very diverse class of the corresponding multi-field interactions in the AT — 1 SYM the¬ 
ories. A rather general class of interesting multi-matrix models with tree-like interactions 
will be studied in the next section. 
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7. Possible generalizations 

A large class of solvable matrix models0 can be classified by ’’tree diagrams”, where 
each edge of the tree connects two matrices sitting at the vertices. The corresponding 
matrix model potential can be written as 


Q 


i/(*) = ■& — J2 +E w ‘(*o 


(7.1) 


i>j=l 


where C V] — 1, if i. j are the neighboring vertices of the tree, and C t] — 0 otherwise. A 
particular kind of such a solvable model with tree-like interaction, the Potts model on 
planar graphs, was first considered in [[42 . 


It is easy to integrate out the “angular” parts of hermitean matrices i = 1,..., Q, 
since they are independent in the case of a tree interaction, taking the corresponding HCIZ 
integrals and to rewrite the partition function of the model in terms of their eigenvalues 
<f>i = diag (z [ l \..., z$) (see for example P5fl ) 


Z = 




k =1 \ i 




i =1 


n A ( ,«) 


2—mi 


(7.2) 


where rn-i — XljLi Cij i s the coordination number of the 4-th vertex. 

Introducing the resolvents of matrices Gi(z) (having as usual, the assymptotics 
Gi(z) — to/ z at z —^ oo) we can write down the following saddle point equations, gener¬ 
alizing the eqs. (p?T5|) of the H2MM: 


Q 


J2 Cij z& - W'(z^) = (m, - 2 )Gi(zW) 

3 = 1 


(7.3) 


As in the case of the two matrix model, this system of equations should be degenerate, 
and this degeneracy is the condition of its solution in terms of an algebraic hyper-surface 
relating all Q variables^. Namely, it should exist a polynomial function of F depending 
on all z^\ i — 1,..., Q, such that 


F (*«...,*«») = 0 


(7.4) 


in analogy with (TT). If the system ( |7.8|) was not degenerate it would give only pointlike 
distributions, leading to the collaps of eigenvalues into one or a few points. 


11 where by “solvability” we mean a possibility to reduce the number of degrees of freedom 
from the order of N 2 to the order of N, by integrating over the angular variables of the matrices 


12 


An important paricular case of such a surface will be considered in [44]. 
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To build the function ( |7.4|) , and to analyse the structure of the corresponding algebraic 
surface we should start as usual, from the “classical” equations (cf. with (|2.2|) ) correspon- 
ing to putting all the r.h.s. of ( |7.3|) to zero. The “classical” limit of the function ( |7.4| ) 
corresponds to the product of all these “classical” equations (in analogy with ( |2.4|) for the 
two matrix model). Then one can write: 


F[z^,...,z^ = JJ 

i= 1 


Q 

E 

3 = 1 


Cij z w - w'(z 


+ deformations = 0, 


where by the deformations we mean adding a polynomial in all variables of lower de¬ 
gree, governed by corresponding multidimensional Newton polyhedron. The coefficients 
in front of the monomials of higher degrees are determined by the assymptotics at infini¬ 
ties following from (|77^) and coincide with their values in the classical part. The rest of 
the deformation coefficients will provide the new moduli of the complex structure of this 
algebraic manifold. 

The algebraic equation ( |7.4| ) has the degree deg — n^Li(-^* — 1), where K, is the 
highest power of the potential W r {z). This corresponds to the number of extrema in 
the classical multi-matrix potential and to the number of moduli parameters of the curve 
contained in the deformation. 

The fact that this algebraic surface is in general not a curve does not contradict the 
existence of the resolvents Gi(z) : which means that for any two variables z^ l \ zS-^ it should 
exist an algebraic curve F^ (z W, z^) — 0 with a polynomial function Fij(x : y ) of the same 
degree. This loss of information in the abovementioned surface might be related to the 
existence of higher dimensional holomorphic differential forms, like the well known three- 
form dil : > on the Calabi-Yau 3-folds. On particular 3-cycles, after the integration over two 
variables it turns into a meromorphic one-differential. This differential is related to the 
resolvent with respect to the third variable [15]. In principle, we could try to restore the 
individual algebraic curves out of the general surface by excluding the variables one by one 
from the equation for the surface, using ( l7..'fj ) . 

Of course it is only a sketch of the construction of the algebraic curve of the tree-like 
multi-matrix model. It would be interesing to precise the details and the structure of these 
algebraic curves, although it might be difficult to do it to the same extent of explicitness 
as we have done in this paper for the two matrix model. 

Some further generalisations are possible if we substitute the potential ( |7.I|) by 


Q 


n«) = T >- - E c , «(* i *j)+E w i(*i) 


(7.5) 


i>j = l 
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where Cij(M ) are arbitrary polynomial functions (nonzero only on a tree). In this case 
we can achieve the reduction to the eigenvalues by the method of character expansion (see 
43(1 and references therein). Namely, we can expand the exponent of each interaction term 


into the GL(N) characters xr{M) 


-C tj Tr _ 




R 


and then use the orthogonality property of the matrix elements to integrate out the relative 
angle of two matrices. An example of such calculation was done in j4(| for the Tr 
interaction. The corresponding system of polynomial equations will include, along with the 
eigenvalue variables, the dual variables corresponding to the highest weights of the Young 
tableaux of the SL(N ) irreducible representations. Here the construction of the algebraic 
surface looks even more difficult, but certainly not impossible. Note that the model still 
stays “solvable” by the character expansion method if we change the arguments of products 
Cij(3>i$j) by with arbitrary integers /c,y , n %3 for each zj-link. 

These models (and those which can be reduced to them by introducing some gaussian 
matrix integrations, like in the case of the Potts model on planar graphs j|2|]) exhaust the 
list of “solvable” (i.e., reducible to the eigenvalues) multi-matrix models. 


8. Conclusion 


In this paper we studied the multi-support solutions of two matrix models and we 
found that they lead to the appearence of a new nontrivial one-dimensional complex ge¬ 
ometry. The corresponding quasiclassical tan-function can be still defined in a standard 
way and even rather explicit formulas defining the tau-funcion can be written down. 

The main problem for the multi-support solutions is nevertheless to write down ex¬ 
plicitly the system of integrable equations which is solved by the corresponding free energy. 
For the one-support solutions such system necessarily includes the dispersionless Hirota 
equations which have a nice and well-known dispersionfull analog. 

The only known analogs of the dispersionless Hirota equations for the multi-support 
case are the associativity or WDVV equations [|47|], and a wide class of their solutions is 
constructed on the base of quasiclassical tau-functions. The formula ( |4. 1 1| ) in particular 
suggests that the quasiclassical partition function of the cubic two matrix model in the 
’’symmetric ansatz”, say with to — 0, satisfies the WDVV equations [47],[48]. This follows 


from simple counting of the number of variables (six coefficients of the potential and three 
filling numbers altogether give the same number of free parameters as the number of critical 
points ( |4.11| ) for n — 2) and the structure of the residue formula for this case. We mention 
this fact since it seems to be the only explicitly known solution to WDVV equations coming 
from the nonhyperelliptic curves (see []40|| for a general discussion of this issue). 
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In this paper we presented a two-support solution of the two matrix model for the 
cubic potential in a rather explicit form. Nevertheless, in principle it can be precised 
further. In particular, the period integrals (|5.37|) are elliptic (according to the structure 
of the underlying curve having the genus g rc< i = 1) and can be in principle calculated 
explicitly. It might be instructive to do it and to write down explicitly their asymptotic 
expansion. 

The multi-solutions we found are quite interesting from the point of view of statistical- 
mechanical models on planar graphs. An example of such model (a double-phase Ising 
model) is described in section 2. It would be interesting to classify all such possible models 
emerging from the multi-support two matrix model. 

From the point of view of the underlying J\f — 1 SYM theory it is very desirable to 
study the degenerations of a higher genus algebraic curve (with more than two cuts filled) 
of the two matrix model considered above and classify the emerging physical excitations 
(monopoles, dyons etc), by analogy to the hyperelliptic solution in (generalized) Seiberg- 
Witten picture. 

Finally, the models described in section 7 should contain a much richer variety of pos¬ 
sible algebraic surfaces describing their planar limit. They certainly deserve a considerable 
attention. The matrix models, due to their natural integrability properties, could give an 
insight into the structure of algebraic surfaces possessing interesting physical applications. 
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